Nonlocal protons and deuterons oppose to disorder: 
A single-crystal neutron diffraction study of 
§ ; KH0.76D0.24CO3 and a theoretical framework 
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C^ ' Abstract. We show with neutron diffraction that a single-crystal of potassium 

hydrogen deuterium carbonate at room temperature, namely KH0.76D0.24CO3, is 

isomorphous with the KHCO3 and KDCO3 derivatives. Protons and deuterons are 

^ ' not disordered particles located at definite sites and Bragg peaks are best fitted 

O . with separate H and D sublattices. We propose a theoretical framework for nonlocal 

observables and macroscopic states compatible with the crystal field. 
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1. Introduction 

Defect-free crystals are macroscopic quantum systems with spatial periodicity. This is 
attested by electronic states and, for nuclei, by discrete phonon spectra observed at any 
temperature below melting or decomposition. However, the quantum nature of nuclei 
is overlooked in many cases. For example, in crystals containing 0-H- ■ -O hydrogen 
bonds, when the coexistence of two configurations for protons, say 01 — H ■ ■ • 02 and 
01 • ■ ■ H — 02, is conceived of as "disorder" of particles with definite positions and 
momenta. This is at variance with vibrational spectra showing that the translational 
invariance of the lattice and the quantum nature of dynamics are not destroyed, even 
on the short time scale of proton dynamics (~ 10~^^ — 10"^*^ s). Such crystals should 
be therefore represented by macroscopic states for nonlocal observables. 

This conflict between disorder and macroscopic states parallels the dichotomy of 
interpretation of quantum mechanics applied to macroscopic systems IH 12]. At the 
microscopic level, a quantum superposition indicates a lack of definiteness of outcome 
evidenced by interferences. Extrapolated to the macroscopic level of everyday life, this 
leads to Schrodinger's Cat in a superposition of \alive) and \dead) states. However, we 
know that a cat is either dead or alive, so we are naturally inclined to think that a 
macroscopic object which has available to it two or more macroscopically distinct states 
is in one of these states for each measurement of the ensemble [31 HI [5]. For hydrogen 
bond protons, this means that dynamics at the microscopic level of a double-well are 
quantum in nature (tunnelling), whereas at the macroscopic level of a crystal, protons 
which have available to them distinct sites are in one of the many possible conflgurations 
for each measurement. 

For open systems, this dichotomy of interpretation can be legitimated by 
decoherence stipulating that an initial superposition state should lose its ability to 
exhibit quantum interferences via interaction with the environment [6J. However, since 
the quantum theory does not predict any dividing line between quantal and classical 
regimes, macroscopic quantum behaviour is possible for systems decoupled from, or 
very weakly coupled to, the surroundings [7]. In principle, there is no upper limit in 
size, complexity, and temperature, beyond which such systems should be doomed to 
classicality. 

The crystal of potassium hydrogen carbonate (KHCO3) is a counterexample to the 
dichotomy of interpretation of quantum mechanics at the macroscopic level. Neutron 
diffraction and vibrational spectroscopy provide positive evidences of macroscopic 
superposition states for protons featuring quantum correlations, at the scale of 
Avogadro's constant, and about room temperature [HI El [IDl [HI [121 [13] . The particle-like 
representation of protons possessing properties in their own right (position, momentum, 
spin) must be abandoned. Nonlocal observables and macroscopic wave functions are 
best appropriate to the translation invariant symmetry of the crystal. These states are 
decoherence-free and there is no transition to the classical regime, as long as the crystal 
structure is not disrupted. 
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The question at issue in the present paper is as to whether or not macroscopic states 
are destroyed by static disorder in a mixed crystal, KHpDi_pC03, as grown from a water 
solution containing H2O, HDO, D2O molecules. In this crystal, H and D nuclei are not 
mobile and there is no segregation into separate domains. Our purpose is to determine 
with neutron diffraction the impact of static disorder to the crystal symmetry. We are 
not aware of such measurements ever published for this system. 

In Section [21 we report neutron diffraction measurements showing that the crystal 
of KH0.76D0.24CO3 is isomorphous with KHCO3 and KDCO3. This is at variance with 
static disorder. Data are best fitted with distinct H and D sublattices. We propose, in 
Section [3], a theoretical framework for superposed proton and deuteron states. 

2. Crystal structures 

2.1. KHCO3 and KDCO3 

Up to room temperature, these monoclinic crystals are isomorphous: space group 
P2i/a (Cg^) with four equivalent KHCO3 or KDCO3 entities per unit cell (Fig. [1]) 
p^ [T5] . They are made up of centrosymmetric diniers (HC0j)2 or (DC0j^)2 linked 
by moderately strong hydrogen (deuterium) bonds, with lengths -Rqho ~ 2.58 A 
(-R0D0 ~ 2.61 A). Dimers lie practically in (103) planes and hydrogen (deuterium) 
bonds are virtually parallel to each other. These crystals are uniquely suited to 
probing dynamics along directions x, y, z, parallel to the stretching (i/OH/z/OD), the 
in-plane bending {60B./60D), and the out-of-plane bending (7OH/7OD) vibrations, 
respectively. 

From 14 K to 300 K, there is no structural phase transition for either crystal. 
The increase of the unit cell dimensions and of -Rqho or Rqbo are marginal, while the 
populations of proton (deuteron) sites change significantly. Below ~ 150 K, only those 
sites corresponding to a unique configuration for dimers (say L, see arrows in Fig. [1]) are 
occupied. At elevated temperatures, less favored R sites, at ~ 0.6 A from the L ones, 
are progressively populated. The space group symmetry and quantum correlations are 
unaffected. This accords with a superposition of macroscopically distinct \L) and \R) 
states corresponding to the P2i/a configurations sketched in Fig. [2] [TT j [T2l[T3] . 

2.2. KHpDi_pC03 

An approximately cubic specimen (3x3x3 mm^) was cut from a crystal obtained by 
cooling slowly a saturated solution in a mixture of H2O (~ 75%) and D2O (~ 25%). 
We utilized the four-circle diffractometer 5C2 based at the Orphee reactor (Laboratoire 
Leon-Brillouin) [16]. The crystal was wrapped in aluminum and maintained at (300 ±1) 
K. The data reduction was carried out with CRYSTALS [l7j. 

Inspection of intensities for absent reflections {OkO ioT k = 2n + 1 and hOl for 
h = 2n + 1) confirms the P2i/a space group assignment. There is no evidence of 
any symmetry breaking due to a static distribution of protons and deuterons. This 



Nonlocal protons and deuterons 4 

information is unique, because neutrons are scattered by nuclei. The same space group 
assignment, should it be determined with X-ray diffraction, would not rule out a static 
distribution of isotopes carrying the same number of electrons. 

We refined a structural model comprising four independent sublattices, H^, H/j, 
Di, D/j, as: 

B.I){ph,Pd,Plh,Prh,Pld,Prd) = q>, 

PhUplh^l + pRH^iR] + PdIpld^l + Prd^r]- 

The sublattices are separate, although they may correspond to identical sites. The 
coherent scattering lengths were set to either bn ~ —3.741 fm for H^, and H/^, or 
bo ~ 6.671 fm for D^, and D/{. Positional and anisotropic thermal parameters for each 
sublattice, isotope concentrations {ph,Pd) and occupancies {plh, Prh, Pld, Prd) were 
allowed to vary independently from each other. 

In Tables [1] to [3l the refined parameters are compared to those of the fully 
hydrogenated analogue at the same temperature. Atoms are labelled H(l), D(l) for L 
and H(2), D(2) for R. The unit cell parameters and the parameters for heavy atom sites 
are virtually identical. Small differences for proton positions are insignificant compared 
to the thermal ellipsoids. 

In the isotope mixture, positional and thermal parameters of protons and deuterons 
at the same sites are identical, so Fig. [1] is representative of the structure with 
a superposition of indiscernible H and D sublattices. In Table [2], the estimated 
concentrations, ph = 0.76(2) and po = 0.24(2) accord with those of the initial water 
solution and the sum piH + Prh + Pld + Prd = 0.99(4) is self-consistent. 

The occupancy prh = 0.19(2) could be slightly smaller than 0.22(1) for KHCO3, 
while Prd = 0.15(2) could be a little bit greater than 0.12(1) for KDCO3 [IS]. A small 
variation of the effective energy difference between L and R configurations with the 
isotope content cannot be excluded, but statistical errors are not conclusive. Finally, 
the O- ■ -O bond length of 2.587(1) A is equal to PhRoko + PdRobo- The lattice of 
heavy atoms could be also made up of separate sublattices, but we did not pursue this 
option. 

For the sake of completeness, we tried an alternative model supposing that H 
and D nuclei merge into hybrid entities [HpDi_p] with the averaged scattering length 
bp = pbu + (1 — J5)6d- These hybrids were incorporated in sublattices [HpDi-p]^ and 
[IIpDi_p]/j as: 

[HpDi_p](pL,pfi) = pi[HpDi_p]L + pR[HpDi_p]K. (2) 

In this model the isotope concentration p is not a free parameter. We set p = 0.76 
and 6o.76 ~ —1.242 fm. When p^ and p^, as well as positional and thermal parameters 
for the L and R sublattices, were allowed to vary independently, the refinement did 
not converge to any stable structure. We imposed further constraints by setting equal 
thermal parameters for L and R. Then, refinement of the remaining 56 free parameters 
did converge to a stable structure (R = 0.036), roughly similar to the previous one. 
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Therefore, this model is not formally excluded by Bragg diffraction. In the next section, 
we present theoretical arguments liable to discriminate models ([1]) and ([2]). 

3. Theory 

The theoretical framework is a follow up of that presented in Refs p!T| [T3] . We begin 
with macroscopic states in pure crystals and then we consider isotope mixtures. 

3.1. The adiabatic separation 

Within the framework of the Born-Oppenheimer approximation, the nuclear 
Hamiltonian can be partitioned as 

Tiv = 'Hn + Hb + Ch,D + 'Hat + CH,D,at, (3) 

where Hu, Hb, and H^t represent the sublattices of H+, D+ and heavy atoms, 
respectively. Ch,d is the coupling term between Tin and Hb, while Cn.o.at couples the 
three subsystems. 

Vibrational spectra show that H and D states are separable, so Ch,d can be ignored. 
The "hybrid" model ([2]) leading to an effective oscillator mass m-np + rri£)(l — p) and a 
continuous frequency shift with p is therefore inappropriate. 

For OHO or ODO bonds, coupling terms between O- ■ -O and H or D degrees of 
freedom are rather large p^lTO], hence beyond the framework of the perturbation theory. 
Yet, quantum dynamics can be rationalized within the framework of the adiabatic 
separation of Tin and TCb, on the one hand, and Hat, on the other [HI |13]. The 
observation of separable H and D states is an evidence of the adiabatic separation from 
heavy atoms. Then, bare protons are fermions while bare deuterons are bosons. 

3.2. Macroscopic states 

Consider a crystal made of N^, A^b, Nc (A/" = Ns,N]^Nc) unit cells labelled j, k, I, along 
crystal axes (a), (6), (c), respectively. The two dimers per unit cell are indexed as j, k, I 
and j', k, I, respectively, with j = j'. For centrosymmetric dimers, there is no permanent 
dipolar interaction, so interdimer coupling terms are negligible [20l [2T] . The eigen states 
of the sublattices of protons or deuterons can be therefore represented with eigen states 
for isolated dimers 

A dimer (H1,H2) or (D1,D2) is modelled with coupled centrosymmetric collinear 
oscillators in three dimensions, along coordinates {aijki} and {a2jki} (« = x,y,z). The 
center of symmetry is at {aojki}- Within the framework of the harmonic approximation, 
the mass-conserving symmetry coordinates independent of j, k, I, and their conjugated 
momenta, 

1 1 

as = —7= (ai - ^2 + 2ao), Psa = —7= {Pla - P2a), 

Y Y (4) 

aa=^(ai + a2), Ps^ = -7=iPla + P2a), 
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lead to uncoupled oscillators at frequencies hjJsa and hw^ai respectively. The difference 
{hjJsa — ^aa) depcuds ou the coupling term (say Aq,). The wave functions, {^njkil'^a)}, 
{^n'jki('^s — "\/2tto)}, cannot be factored into wave functions for individual particles, 
so there is no local information available for these entangled oscillators. Then, the 
wave functions for protons or deuterons are subjected to the symmetrization postulate 
of quantum mechanics in different ways [22]. As the wave functions {^njki('^a)} and 
{^n'jki('^s — v^tto)} are unaffected by permutation changing the signs of {aojki}, they 
are appropriate for bosons, but not for fermions. 

3.2.1. Protons In order to antisymmetrize the degenerate ground state with respect 
to permutation, the wave functions are rewritten as linear combinations 

©ojkii = -7^\{ ^ojki(«a) [^ojki(«s - V2ao) ± ^ojki(«s + V2«o)] (5) 

and the state vectors with singlet-like IS") or triplet-like \T) spin symmetry are: 

|eojki+)®|5); ,s 

|eojki-)®|T). ^^ 

The oscillators are now entangled in position, momentum, and spin. As there is no level 
splitting, entanglement is energy-free and independent of A^. 

The spatial periodicity of the sublattice leads to collective dynamics and nonlocal 
observables in three dimensions. With the vibrational wave function for the unit cell 
J, fc, /, namely Hojkir = ©ojkir + tQoj'Mt, where r = "+" or "— " for singlet-like or 
triplet-like symmetry, respectively, we write phonon waves as 

. Nc Nb Na 

Sor (k) = -j= ^^Y^ Sojkir exp (ik ■ L) , (7) 



^ 1=1 k=l j=l 

where k is the wave vector, L = ja + A;b + /c is the lattice vector, a, b, c, are the unit 
cell vectors. Then, antisymmetrization leads to 

k ■ L = modulo 2ti. (8) 

This means that there is no phonon (no elastic distortion) in the ground state ("super- 
rigidity") [11]. The lattice state vectors 

|H+) = |So+(k = 0))®|5) 

|H_) = |5o-(k = 0))®|T) ^ ' 

are superposed in the ground state as: 

v^|So+(k = 0))®|^); 

VAr|So-(k = o))®|r). ^ ^ 

These states are intrinsically steady against decoherence. The main source of 
disentanglement is the thermal bath but, even at room temperature, the population 
of the first excited state (< 1% for 7OH ^ 1000 cm~^) is negligible. 
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3.2.2. Deuterons For bosons, the wave functions for a dimer 

©Ojkl = n ^Ojkl(«a)^ojkl(«s - V^«0), (11) 

a 

are symmetric with respect to particle permutation. As there is no significant spin-spin 
couphng compared to the thermal bath, the spin of a pair is statistically zero, so there 
is no correlation at all. Phonon states can be written as 

. N, Nb Na 

l^i^i) = ^17 E E E l-0Jki±) exp(zk ■ L), (12) 

vA/ 1^1 k=i j=i 

where Sojki± = ©ojki ± ©oj'ki- In contrast to the fermion case, there is no super-rigidity. 
Needless to say, the H and D sublattices have the same number of degrees of freedom (12 
per unit cell) , but the symmetrization postulate shrinks the size of the allowed Hilbert 
space from ~ 12 for bosons to ~ 12A/' for fermions. 

3.2.3. Isotope mixtures The sublattice at room temperature can be represented 
with a superposition of separable proton and deuteron states, each of them being a 
superposition of L and R states, as: 

pf[{pTH\^±)L + P«/J|H±)«] + p'J^[pTD\^^.L±)L + pil\D^R±)Rl (13) 

3.3. Neutron diffraction 

Neutrons (spin 1/2) are unique to probing the spin-symmetry of proton states (j9]). 
However, the spin entanglement is extremely fragile, so only "noninvasive" experiments, 
free of measurement-induced decoherence, are appropriate |1] . This means: (i) no energy 
transfer; (ii) no spin-flip; (iii) particular values of the neutron momentum transfer vector 
Q [23] preserving the super-rigidity. This is realized when the components Qx, Qy, Qz, 
match a node of the reciprocal sublattice of protons pT| [13] . 

With a reactor-based diffractometer, these conditions are not realized at the 
measured Bragg peaks. The spin-symmetry is destroyed, the super-rigidity is lost, and 
the final proton states |Hk±) are analogous to boson states (fT2l) . The structural model 
([1]) is therefore consistent with diffraction by an initial superposition state flT3l) and one 
of the final states |Hk-|-)i, |Hk±)i?, |Dk±)L, |Dk±)ij. The final state of a particular event 
is unknown because the four states are degenerate in reciprocal space. The probabilities 
of these scattering events are given by ([1]), in accordance with (fT3l) . 

4. Conclusion 

The space group symmetry of KH0.76D0.24CO3 determined by single-crystal neutron 
diffraction shows that H and D nuclei are not individual and mutually exclusive particles 
located at definite sites. There is no disorder, either static or dynamic. Bragg peaks are 
best fitted with separate H and D sublattices of identical sites. Each sublattice is further 
made up of distinct L and R sublattices. This structure accords with a superposition 
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of separable macroscopic states. Further measurements should provide information as 
to whether the spin-symmetry of proton states survives in isotope mixtures. 

There is now a substantial amount of experimental data showing that crystals 
of KHCO3 and isotope derivatives are macroscopic quantum objects on the scale of 
Avogadro's constant and room temperature. The wave, or matter-field, representation 
is imposed by the electronic structure insensitive to isotope substitution. The theoretical 
framework is based on fundamental laws of quantum mechanics. It is free of any ad 
hoc hypothesis or parameter. There is every reason to suppose that macroscopic states 
should occur in many systems. Neutron diffraction by partially deuterated hydrogen 
bonded crystals is well-suited to evidencing the nonlocal character of protons and 
deuterons, as opposed to disorder. 
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Table 1. Neutron single crystal diffraction data and structure refinement at 300 K. A 
= 0.8305 A, space group P2i/a. Tire variance for the last digit is given in parentheses. 



Crystal data 


KH0.76D0.24CO3 


KHCO3 [llj 


a(A) 


15.180(1) 


15.180(1) 


6(A) 


5.620(1) 


5.620(4) 


c(A) 


3.710(1) 


3.710(4) 


P 


104.67(1)° 


104.67(5)° 


Volume 


306.2 


306.2 


Reflections measured 


1288 


1731 


Independent reflections 


1011 


1475 


Reflections used 


811 


1068 


(t(I) limit 


3.00 


3.00 


Refinement on F 






R-factor 


0.026 


0.035 


Weighted R-factor 


0.038 


0.034 


Number of parameters 


77 


56 


Goodness of fit 


1.043 


1.025 


Extinction 


121.8(40) 


3260(100) 
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Table 2. Atomic positions, isotropic temperature factors and site occupancies for 
KH0.76D0.24CO3 (first lines) and KHCO3 (second lines) [TT] at 300 K. The variance 
for the last digit is given in parentheses. 



Atom 


x/a 


y/h 


z/c 


U(iso)(A2) 


Occupancy 


K(l) 


0.16534(6) 


0.0242(2) 


0.2956(2) 


0.0232 


1.000 




0.16534(6) 


0.0230(2) 


0.2952(2) 


0.0213 


1.000 


C(l) 


0.11963(3) 


0.51735(8) 


-0.1443(1) 


0.0185 


1.000 




0.11957(3) 


0.51636(8) 


-0.1444(1) 


0.0164 


1.000 


0(1) 


0.19345(4) 


0.5321(1) 


0.0947(2) 


0.0279 


1.000 




0.19344(4) 


0.5307(1) 


0.0946(2) 


0.0259 


1.000 


0(2) 


0.08237(4) 


0.3220(1) 


-0.2734(2) 


0.0260 


1.000 




0.08237(4) 


0.3206(1) 


-0.2735(2) 


0.0238 


1.000 


0(3) 


0.07737(4) 


0.7195(1) 


-0.2741(2) 


0.0266 


1.000 




0.07745(4) 


0.7186(1) 


-0.2743(2) 


0.0243 


1.000 


H(l) 


0.0158(3) 


0.6897(3) 


-0.4492(3) 


0.0310 


0.643(7) 




0.01631(12) 


0.6905(2) 


-0.4491(5) 


0.0330 


0.823(4) 


D(l) 


0.0158(3) 


0.6897(3) 


-0.4492(3) 


0.0310 


0.200(6) 


H(2) 


-0.019(2) 


0.686(2) 


-0.560(2) 


0.0324 


0.12(1) 




-0.0207(6) 


0.680(1) 


-0.563(2) 


0.0338 


0.177(4) 


D(2) 


-0.019(2) 


0.686(2) 


-0.560(2) 


0.0324 


0.03(1) 
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Table 3. Thermal parameters in A^ units for KH0.76D0.24CO3 (first lines) 
and KHCO3 (second lines) [H] at 300 K. The variance for the last digit 
is given in parentheses. These parameters account for the contribution 

of each atom to Bragg's peak intensities through the thermal factor T"* = 
exp [-27r2(C/f/h^a*2 + t/fik^^*^ + ^fll^c*^ + 2[/f2*ha*k6* + 2U^l'kh*\c* + C/fi*lc*ha*)] , 
where a*, b*, c*, are the reciprocal lattice parameters and h, k, 1, the indexes in 
reciprocal space. 



Atom 


Un 


U22 


U33 


U23 


Ui3 


U12 


K(l) 


0.0265(4) 


0.0231(4) 


0.0181(3) 


0.0001(3) 


0.0018(3) 


0.0006(3) 




0.0240(3) 


0.0222(3) 


0.0159(3) 


-0.0004(3) 


0.0017(2) 


0.0009(3) 


C(l) 


0.0178(2) 


0.0194(2) 


0.0171(2) 


-0.0003(1) 


0.00207(13) 


-0.0003(1) 




0.0150(2) 


0.0182(2) 


0.0147(2) 


-0.0000(1) 


0.0016(2) 
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Figure 1. Schematic view of the crystahine structure of KHCO3 at 300 K. The arrows 
point to the sites occupied at low temperature (the L configuration, see text). Dashed 
Hues through protons are guides for the eyes. The eUipsoids represent 50% of the 
probabiUty density for nuclei. 
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Figure 2. Projections on the (a, b) plane of the L (bottom) and R (top) configurations. 



